Clinics and Practice 2013; volume 3:e30 



Correspondence: Shigeo Tanabe and Yoshihiro 
Muraoka, 98-1 Dengakugakubo, Kutsukake, 
Toyoake, Aichi, 470-1192 Japan. 
Tel. +81.562.93-9000 - Fax: +81.562-95.6817. 
E-mail: tanabes@fujita-hu.ac.jp 

Keywords: gait analysis, gait trajectory, treadmill. 

Contributions: ST, YM, YT, study execution, data 
analysis, manuscript writing; ES, KO, TT, study 
execution and design; all authors discussed 
results and commented on the manuscript. 

Conflict of interests: the authors declare no 
potential conflict of interests. 

Received for publication: 10 July 2013. 
Revision received: 21 September 2013. 
Accepted for publication: 24 September 2013. 

This work is licensed under a Creative Commons 
Attribution Noncommercial 3.0 License (CC BY- 
NC 3.0). 

©Copyright S. Tanabe et ai, 2013 
Licensee PAGEPress, Italy 
Clinics and Practice 2013; 3:e30 
doi:10.4081/cp.2013.e30 



Simple method to reduce 
the effect of patient 
positioning variation on 
three-dimensional motion 
analysis during treadmill gait 

Shigeo Tanabe,^ Eiichi Saitoh,^ 
Kei Ohtsuka,^ Toshio Teranishi,^ 
Yutaka Tomita,^ Yoshihiro Muraoka^ 

^Faculty of Rehabilitation, School of 
Health Sciences, Fujita Health University; 
^Department of Rehabilitation Medicine I, 
School of Medicine, Fujita Health 
University; ^Fujita Memorial Nanakuri 
Institute, Fujita Health University, Japan 



Abstract 

Recently, three-dimensional (3D) closed 
curve trajectories of markers placed at strate- 
gic body locations, called cyclograms or 
Lissajous-like graphs, are used for treadmill 
gait analysis. A simple method is presented to 
reduce the effect of patient positioning varia- 
tion. After breaking down movement into three 
components (anterior-posterior, medial-lateral 
and superior-inferior), the time-series data 
and time-inverted data are serially concatenat- 
ed. A fast Fourier transform (FFT) is done, and 
a high-pass filter (except 0 Hz) is applied to 
the anterior-posterior and medial-lateral com- 
ponents. Next an inverse FFT is executed, and 
the posterior half of the outcome, correspon- 
ding to time-inverted data, is deleted. The 3D 
closed curve is then reconstructed. Results 
showed that the proposed method was able to 
reduce the effect of patient positioning varia- 
tion. Since the adjusted curve is simply a sym- 
bolized gait pattern, the method might be use- 
ful as an adjunct tool in observational gait 
analysis. 



Introduction 

Bipedal ambulation is one of the most com- 
mon human physical activities. In the clinical 
setting, many patients with neurological or 
orthopedic disorders demonstrate distur- 
bances in gait. Gait exercise on a treadmill is 
effective rehabilitation after neurological 
assault.^ ^ In addition, treadmill exercise has 
been combined with robotic assistance and/or 
functional electrical stimulation."^ In order to 
evaluate the effects of exercise and to deter- 
mine optimal exercise conditions, it is neces- 
sary to objectively clarify the ambulation 
capacity of each patient. Some spatiotemporal 



parameters (swing time, stance time, double 
support time, etc.) measured by a motion 
analysis system have been used in previous 
studies.^^ Recently, three-dimensional (3D) 
closed curves of the trajectories of markers 
placed on body landmarks, called cyclograms or 
Lissajous-like graphs, have been used for gait 
analysis.^^^ Since the measurements are 
acquired on a treadmill, each marker returns 
to almost the same position after one stride 
cycle, and consequently, the path forms a 
closed curve. In general, the curves are used 
after decomposing them into three planes 
(frontal, sagittal, and horizontal). Because the 
curves indicate human motion patterns, physi- 
cians and therapists can interpret gait patterns 
from the curves. However, it is difficult to apply 
this study method to gait disorders because 
this methodology assumes that patients walk 
in a consistent position which cannot be guar- 
anteed. Even in healthy persons, the position- 
ing of anterior-posterior and medial-lateral 
directions vary from stride to stride, thus caus- 
ing drifts of curves (Figure 1).^^ 

The drifts can distort the curve shape and 
increase the standard deviation of averaged 
curves. Such drifts might obscure the specific 
or pathognomonic points of curves. 
This paper presents a simple method to reduce 
the effect of patient positioning variation on 
3D motion analysis data during treadmill gait. 
The method uses a fast Fourier transform 
(FFT) and high-pass filter, and is expected to 
help motion analysis using 3D closed curves in 
gait disorders. 



Materials and Methods 

Ethics approval 

The ethics committee of the Fujita Health 
University approved the study protocol (No. 09- 
028). The experiment was conducted in accor- 
dance with the 1975 Declaration of Helsinki 
(rev 2008). 

Experimental procedure 

One healthy young subject with no history of 
neurological or orthopedic disorders, who was 
expert with substantial experience treating 
various gait disorders, participated in the 
experiment. Prior to data collection, the sub- 
ject gave written informed consent after thor- 
ough explanation of the study. The subject 
walked on a treadmill (Adal 3D, Medical 
Development, Nantes, France) with a colored 
marker placed on the right lateral malleolus. 
The subject was instructed to mimic a gait dis- 
order defined as fluctuation in the positioning 
in the anterior-posterior and medial-lateral 
directions. The belt speed was set at 4.0 km/h. 
Before obtaining measurements, practice tri- 



als were allowed for the subject to become 
familiar with a disordered gait. After familiar- 
ization, the marker trajectory was collected by 
a motion capture system (Motion recorder, 
Kissei Comtec Co., Ltd, Nagano, Japan) at a 
sampling frequency of 60 Hz. Standard data 
processing was performed using KineAnalyzer 
software accompanying the device. 

Processing method 

The proposed method is composed of five 
steps. First, after decomposing the motion into 
three components (anterior-posterior, medial- 
lateral and superior-inferior), each time-series 
data and time-inverted data are serially con- 
catenated because the method uses a FFT algo- 
rithm that assumes that the signal is periodic, 
thus beginning and ending with the same 
value, otherwise spectral leakage artifacts 
would occur. Second, FFT computation with a 
rectangular window function is used to obtain 
the frequency spectrum. Third, a high-pass fil- 
ter is applied to the anterior-posterior and 
medial-lateral components of the data (not 
applied to the superior-inferior component). 
Specifically, part of the low frequency band [up 
to one-third of the gait cycle frequency (except 
0 Hz) ] is replaced with zeroes since this where 
the positioning fluctuation occurs. Fourth, the 
inverse FFT computation reconverts the fre- 
quency spectrum to time-series data. Finally, 
after deleting the posterior half of the data 
(corresponding to inverted time-series data), 
the 3D closed curve is reconstructed. 



[page 84] 



[Clinics and Practice 2013; 3:e30] 



OPEN /jl ACCESS 



Brief Report 



Results 

Figure 2 shows the waveform at each pro- 
cessing step in the sagittal plane. 

In the first step, the first and last points of 
the waveform became the same value by cou- 
pling to the inverted time-series data. In the 
third step, the proposed method selectively 
removed noise content (low frequency band) 
resulting from positioning fluctuation. 

Figure 3 shows the closed curves in the 
sagittal and the horizontal planes, respectively. 

The noise caused by fluctuation in position- 
ing in the anterior-posterior and medial-lateral 
directions was removed by application of the 
proposed method. 



Discussion 

The results suggest that the proposed 
method is able to reduce the effect of patient 
positioning variation on the 3D closed curve 
measured during treadmill gait. In theory, vari- 
ations of each marker made up of two inde- 
pendent components. One is the variation of 
gait pattern itself. It is expressed in per-cycle 
variation of marker position relative to center 
of gravity. The other is the variation of patient 
position. It is caused by the variation of gait 
speed (anterior-posterior component) and a 
decrease of straight-walking stability (medial- 
lateral component). Although both components 
have important clinical implications, first step 
in the gait analysis should be to evaluate the 
variation of gait pattern after removal of the 



effect of patient positioning variation. Since 
the two components differ in frequency, the 
latter may be selectively filtered. In particular, 
the frequency of gait pattern variation is high- 
er than gait cycle whereas sufficiently lower in 
patient positioning variation. In the proposed 
method, the gait pattern is cleared up by high- 
pass filter [one-third of the gait cycle frequen- 
cy (except 0 Hz)] for reducing the effect of 
patient positioning variation. 

Previous studies examined the time vari- 
ability of gait parameter in various gait disor- 
ders. Patients with Parkinson's disease have 
greater variability in step time^^ and stride 
time.^"^ Post-stroke patients have variability in 
stride time compared to control subjects.^^ 
Patients with cerebral palsy also display 
increased variability of step time and corre- 
sponding higher center of pressure dispersion 




Figure 1. Typical example of the effect of fluctuation in the posi- 
tioning of anterior-posterior direction on the marker trajectory 
of the right lateral malleolus in the sagittal plane. 



Figure 3. Raw and compensated marker trajectories attached to 
the right lateral malleolus. A) View from the sagittal plane. B) 
View from the horizontal plane. 
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relative to unaffected subjects.^^ Since time 
variability might cause the patient positioning 
variation, the proposed method is especially 
important in measuring gait in those affected 
with these disorders. 

Compensating 3D closed curves with the 
proposed method might hold great promise for 
improving the efficacy of observation gait 
analysis that is widely used in clinical practice. 
Physicians and therapists generally perform 
pattern recognition in gait disturbances to 
plan patients' treatments. However, observa- 
tion analysis does require a higher skill level. 
In detecting gait anomalies in unilateral 
below-knee amputees using observation ver- 
sus biomechanical analysis, the pick-up rate by 
observation was 22.2% when compared with 
biomechanical analysis.^^ When comparing 
visual and 3D gait analysis in patients with 
cerebral palsy, only two out of ten anomaly 
points could be similarly evaluated by both 
methods.^^ Brunnekreef et al}^ suggest that a 
great deal of experience in observation gait 
analysis is required to improve observation 
skills. In other words, some experience is not 
sufficient; only expert raters accomplish sig- 
nificantly better reliability levels of visual gait 
observation compared to inexperienced raters. 
Although the unmanipulated 3D closed curve 
is distorted by the effect of patient positioning 
variation, the compensation provided by the 
proposed method reduces the shape complexi- 
ty to the symbolized gait pattern that has great 
understandability. Therefore, it might be use- 
ful as an adjunct tool in gait pattern recogni- 
tion. 

Some limitations have to be addressed in 
future work. First, large number of repeated 
measurements is required to clarify intra- and 
inter-subject variability. Second, a wide variety 
of disorders should be measured since gait 
abnormality might be affected by type of disor- 
der. Third, the processing method should be 
analyzed in detail. The serially concatenated 
data (time-series data and inverted data) was 
used to conform the last point of the signal to 
the first point. However, the data introduces a 



discontinuity in the first time-derivative at the 
concatenation point. It might have effects on 
an amplitude spectrum of the signal using FFT. 
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